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The magnetic characteristics for superconducting samples of different ge- 
ometries (plates and films of various widths) and orientation with respect to 
external field are compared. A similarity is established between the magneti- 
zation curves M{H) and the corresponding susceptibilities. A sign reversal is 
predicted for the 3rd harmonic of the ac susceptibillity when the irreversibil- 
ity mechanism changes from bulk pinning to edge barrier. A mutual influ- 
ence of the edge barrier and bulk pinning on the magnetization curves and 
field-dependent ac susceptibility is studied on the basis of an exactly solvable 
model (narrow plate in a parallel field and/or narrow film in a perpendicular 
magnetic field) . 



I. INTRODUCTION 

Study of the magnetic characteristics of type-II superconductors and their behaviour in 
a magnetic field varying in time is of active interest from both the fundamental and applied 
standpoints. A large number of theoretical and experimental works [TTJ have been devoted 



to investigation into the magnetization of samples of various geometries, hysteresis losses 
and ac susceptibilities Xn C = Xn + ^Xn- Two mechanisms are known to be responsible for 
these characteristics, specifically, bulk pinning and edge barrier that prevents vortex entry 
into (exit from) a smooth-surface sample. The effects produced by these mechanisms on 



the shape of the magnetization - and/or hysteresis curves have been considered by different 
authors P,f4|,|3|-|TT|| . However, the influence of an edge barrier on field dependences of the 
hysteresis losses and the behaviour of higher harmonics of susceptibility have not been 
studied adequately so far. In particular, as reported in our recent publication |I2"| , the 



allowance for the edge barrier for a narrow film in a transverse magnetic field results in 
a sign reversal of the 3d harmonics of susceptibility ^{H) and x 3 '(ff). This is in stark 
contrast with the situation in which bulk pinning is the major mechanism of irreversibility. 
Therefore, it would be extremely important to analyse a joint influence of an edge barrier 
and bulk pinning on the irreversible characteristics of type-II superconducting samples. For 
bulk superconductors a similar problem was studied in by Clem who has formulated 
a generalized local model of the critical state. This model provides an adequate basis 
for calculation of the magnetization curves and ac susceptibility for bulk superconductors 
13|,[14|. Yet, in low-dimensional superconductors (single crystals with high demagnetization 



factor, films) the effect of edge barrier and bulk pinning on magnetic characteristics is still 
not fully understood. For example, in ]T3| the authors take account of the joint influence of 



edge barrier and bulk pinning on the shape of the magnetization curve for a wide film in a 
perpendicular magnetic field. Unfortunately, the dependences M(H) obtained therein cannot 
be expressed in an analytical form, which does not allow to find the field and temperature 
dependences of magnetic susceptibility. It is, therefore, highly desirable to formulate an 
analytically solvable generalized model of the critical state for superconductors of arbitrary 
geometry, both longitudinal (as shown in Fig. la) and transverse (see Fig. lb) with respect 
to the external magnetic field. 

In this work we report a study on the behavior of magnetization and the harmonics of ac 
susceptibility as a function of the amplitude of an applied low-frequency magnetic field. As 
an example, we used an exactly solvable model of a narrow film with an edge barrier. This 
model accounts for arbitrary bulk pinning which is characterized by a certain dependence of 
the depinning current density j p on magnetic field if. For specific calculations, two model 
dependences j p (H) were used. One of them assumes that the depinning current density j p 



is a constant value independent of H (the Bean model). In the other case the magnetic 
characteristics were obtained for the Kim- Anderson (KA) model dependence j p (H). 

The paper is organized as follows. In Section II we formulate a generalized model of 
the critical state in a narrow film/plate. On the basis of this model we obtain current 
densities and vortex distributions in the mixed state taking into account both edge barrier 
and bulk pinning as major irreversibility mechanisms in type-II superconductors. Two 
specific depinning models were studied: the Bean model (see II. A) and the Kim-Anderson 
model (see II. B). The magnetization curves for a narrow film/plate are obtained and the 
harmonics of ac susceptibility are calculated in Sec. III. In Sec. IV a comparison is made 
of the magnetic characteristics for superconducting samples of different geometries (plates 
and films of various widths) and orientation with respect to external field. A similarity is 
found out between the magnetization curves M(H) and the corresponding susceptibilities. A 
sign reversal effect is predicted for the 3d harmonic of ac susceptibillity, following a change 
of irreversibility mechanisms (from bulk pinning to edge barrier). In Sec. V we discuss 
the applicability conditions of the approach used. The conditions of the barrier-to-pinning 
crossover in the temperature-dependent ac susceptibilities are analysed. A brief summary 
of our results is presented in Sec. VI. 



II. GENERALIZED MODEL OF THE CRITICAL STATE 

Let us consider a narrow plate (Fig. la) or a film (Fig. lb) of width W and thickness 
d < A (W < A for the plate, W < \ 2 /d for the film; A being the London penetration depth), 
which is infinite in the x-direction, placed in a magnetic field H = (0,0, H). In view of 
the symmetry of the problem, the screening current has one component j x = j(y), which 
satisfies the equation (see Appendix) 

Here C = cW/ (8tt\ 2 ), n is the vortex density, O is the magnetic flux quantum, and the 
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dimensionless coordinate \y\ < 1 is measured in units W/2; the origin of the coordinates is 
assumed to be in the sample centre (y — 0). Equation (1) can be derived from the full version 
of the Maxwell-London equation (see e.g. ref. ||) by neglecting the nonlocal (integral) term. 

Since in the absence of transport current n(y) = n(—y), and j(y) = —j(—y), we may 
reduce our consideration to only one half of the film (for definiteness, the right-hand one, 
y > 0). It is assumed that a superconductor is placed in a magnetic field which is first 
quasistatically increasing to some value H and then decreasing down to — H Q . As is known, 
given an edge barrier the vortices penetrate a sample when the near-the-edge current den- 
sity reaches some critical value j s (for ideal surface j s coincides with the Ginzburg-Landau 
depairing current density jcLj- Therefore, up to the field H = H s = j s /C the film remains 
in the Meissner state (see Fig. 2a): 

j(y) = yHC , < y < 1, (2) 

Further increase of the magnetic field will cause penetration of vortices that will move 
inside a superconductor until the Lorentz force F = j(fio/c acting on each vortex becomes 
equal to the pinning force F p = j P 4>o/c. 

A. Field-independent bulk pinning 

First consider a model in which the depinning current density is independent of H (the 
Bean model). The vortices that have penetrated into the superconductor at H > H s will be 
distributed with the density n(y) = H/(p in a region a(H) < y < b(H), where a(H) = H p /H 
defines the depth of vortex penetration, H p = j p /Co. Unlike the classical Bean model, the 
field of complete penetration (a = 0) here is H = oo, which is a general property of the 
nonlocal model of the critical state both in bulk [|T6| , p!7| and quasi-two-dimensional 0,|9| 
superconductors. The parameter b(H) = 1 — (H s — H p )/H determines the size of the vortex- 
free region (of width 1 — b = (H s — H)/H), from which the vortices are "driven" by a strong 
near-edge current j > j p . The corresponding distribution of the current density has the form 
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(see Fig. 2b): 



yHC 0<y<a, 
j(y) = { 3 P a<y<b, (3) 

HC (y-l)+j s b<y<\. 
We now start reducing the magnetic field from H to — H . Until the current density 
in the vortex-occupied region a(Ho) < y < b(Ho) has reached the critical value, j = —j p , 
the distribution of vortices will be "frozen" and the distribution of current density will be 
described by the expression 



yHC < y < ao, 

j(y) = { y(H- H )C + a H C a < y < b , (4) 
yHC + H C {a - b ) b < y < 1, 
where a = a(H ), b = b(H ). The corresponding distributions j(y), n(y) are shown in 
Fig. 2c. 

At the flux- defreezing field H d f = H (H — H s — H P )/(H — H s + H p ), the current 
density j(bo) becomes equal to — j p and the vortices start moving towards the film edge. 
The magnetic flux related to the vortices will stay in the film until the vortex-exit field, 
H ex = H — H s — H p) is reached. So, in the field interval H ex < H < H d f the current density 
and vortex distribution will have the following form (see Fig. 2d) 



j(y) 



yHC < y < a , 

y(H - H )C + a H C a < y < h, 
-j P h < y < b 2 , 

yHC + H C (a -b ) b 2 < y < 1; 



(5a) 
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n{y) = < 



(56) 



< y < a , 
H /% a <y < 6i, 
H/% h<y < 6 2 , 
6 2 < y < 1, 

where 6i = 2H P /(H — H), b 2 = (H — H s — H p )/H. A further decrease in the magnetic 
field starting from H ex will cause the vortices exit from the film at < H < H ex , as well 
(see Fig. 2e): 



j(y) 



yHC < y < a , 

y(H - H )C + a H C a <y <h, 
-j P h < y < 1; 



(6a) 



< y < a, 
n(y) = { H /% a<y <b u (66) 
H/$o 6i < y < 1. 

Note that at H = the remaning vortices occupy only the region a < y < 6 X (see Fig. 
2f). The edge current density in this case exceeds — j s (which prevents entry of antivortices 
into the film), therefore, with a further decrease of the field the distribution 



j(y) 



yHC < y < a , 

y(H - H q )Cq + a H C a <y <h, 
yHC - H C (a - 6i) h < y < 1; 



(7a) 



< y < a , 
n(y) = \ Ho/^o a <y <b 1} 
6i < y < 1, 
will be valid right down to the field = (H 



(76) 



H p )/2 



Hq — H s — Hp) 2 /A + Hq(H s — H p ), at which vortices of the opposite sign ('antivortices') 
will start entering into the film. In the field range —Hq < H < H^ (see Fig. 2g) one finds: 



j(y) 



yHC < y < a , 

y(H - H )C + a H C a < y < b u 

-j P h < y < 6 3 , 

yHC -(H + H s ) 63 < y < i; 



n(y) 



< y < oq, 

H /$o a <y <b u 

H/$o h < y < 63, 

63 < y < 1. 



where fe 3 = 1 + (if a - H p )/H. 



(8a) 



(86) 



B. Field-dependent pinning (Kim- Anderson model) 

Let us find the distribution of vortices and a current in a film with a magnetic-field- 
dependent bulk density of the depinning current j p (H) = C H% 2 /(Hki + \H\) (KA model), 
taking edge barrier into account. We shall follow the same procedure as considered above. 
Thus, for the magnetic field increased to the magnitude H , j(y) is determined by expressions 
(3), where a(H) = H 2 k2 /(H(H kl + \H\)), b(H) = a(H) + (H - H s )/H. When the field is 
decreasing from H to H d f, the current density and vortex distributions are similar to those 
in (4) with new values for a = a(H ), b = b(H ), and the field Hdf being determined from 
the equation 

b (H - H ) + a H + H " 2 = 

tiki + \ti I 

Then, with the field decreasing from H d f to H ex , j(y) and n(y) are defined by the 
dependence (5) in which now we have 

^_ Hki 1 | 1 b ^_ H b + b 1 (H-H ) 

Hq — H Hki + -f^o Hti + \H\ H 
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and field H ex is determined from the condition b 2 (H ex ) = 1. If the field amplitude H 
satisfies the condition H < y/2Hki, then with a further decrease of the magnetic field down 
to — H the evolution of j(y), n(y) distributions is described by formulae (6-8) in which 6 3 
now is expressed as 



H s + H 



H k2 



H H(H kl + \H\Y 
and field satisfies the equation 

H + H (a - h(H)) + H s = 0. 

For the amplitudes that meet the condition Hq > y/2Hki, as can be easily shown, the 
dependence b\(H) posesses a minimum at some field h\ + ^ = (\/2 — 1)(H Q + \/2H k i) > 0. 
In this case expressions (6) will hold only at H > h\ + \ A detailed examination shows that 
with a still further decrease of the field the distribution of the current density and vortices 
density in a film takes the form 



o(y) 



yHC < y < a , 

y(H - H )C + a H C a < y < b*, 

My) b* <y <b 2a , 

-3 P {H) b 2a <y<\] 



n(y) 



(9a) 



< y < oq, 

H /<5>o a <y <h, 

(H kl + H 2a )(b* - y)/2b*% + H 2a /$ b* < y < b 2a , 
H/% b 2a <y< 1, 

where 6* is determined from the condition 6* = bi(H^), b 2a from the equality jo{b 2a ) 
-jp(H), and Mv) ls defined as 



(96) 



Mv) = C (H - Hi +) )y - Co 



tt , rr(+) 
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This current distribution is valid in the field range < H < H[ +) . If the field keeps 
decreasing, then, by analogy with the Bean model (provided H s > H^ 2 /Hki), antivortices 
at first will not be able to enter the film, and the current /vortex distribution in the film will 
be described by the following expressions 



j(y) 



yHC < y < oq, 

y(H - H )C + a H C a < y < b*, 
jo(y) b* <y <b 2a , 

yHC - j p (H) - b 2a HC b 2a < y < 1; 



(10a) 



n(y) 



(106) 



< y < oq, 

H /<S>o a <y < h, 

(H kl + H 2a )(b* - y)/2b*% + H 2a /$ b* < y < b 2a , 

o b 2a < y < l. 

With a further decrease of the magnetic field two situations are possible, depending on 
the problem parameters. In one case, at first b 2a reduces to 6* (which occurs at H — h[ - 1 ), 
following which the current density at the edge reaches the antivortices-entry threshold —j s . 
In the other case, the edge current density first falls down to — j s (it takes place at some 
field H — H), and then b 2a decreases to 6*. For definiteness, we shall consider only the first 
situation, < \H\. Here expressions (10) will hold right down to the field H = h[~^ 

which is determined by the condition b 2a = b*. Further, down to the field H = —H , j(y) 
and n(y) are found from (7,8) with the 6 3 and that have been introduced above. 

Note that a similar method can be used to obtain the current /vortex density distribution 
for an arbitrary function j p (H). For a monotonically decreasing function j p (H) the resulting 
distributions will have a qualitative similarity to those obtained above. In the peak-effect 
conditions characterized by a nonmonotonic behaviour of j p (H) the issue on the form of j(y) 
and n(y) distributions needs to be addressed separately. 
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III. MAGNETIZATION CURVES AND AC-SUSCEPTIBILITY HARMONICS 

FOR NARROW FILMS 



A. Magnetization curves 

The above expressions for the current density allow to find the magnetization dependence 
for superconducting plates (films) in different ranges of field and obtain a hysteresis curve. 
For the Bean model given H > H s + H p the magnetization curve M(H) is defined as follows 



H (b 3 - ajj)/2 + 3H (a - b )/2 + H, 



Hdf < H < H 



M(H) = 7 < 



bl(H - H)/2 - H a 3 /2 + Hb\/2 + H + 3H (a - b )/2, H ex <H< H df , 

b 3 {H - H)/2 - H al/2 + 3b x {H - H )/2 + 3H a /2, 0<H< H ex , 

H {bl - al)/2 + 3H (a - h)/2 + H, < H < 0, 

b\{H - H)/2 - H a 3 /2 + H - 3{H + H s )/2 + Hb 3 3 /2, -H < H < H<f\ 

(11) 

where 7 = -C W/12c. 

At a relatively low field magnitude, H < H s + H p , the mode described by expressions 
(4) will immediately change for the mode (7). Thus, the dependence M(H) in the fields 
interval H s < Hq < H s + H p will have the form 



M(H) = 7 < 



H (bl - al)/2 + 3H (a - b )/2 + H, H df < H < H , 

H Q (bl - al)/2 + 3H (a - h)/2 + H, <H< H df , 

b\{H - H)/2 - H al/2 + H - 3{H + H s )/2 + Hb 3 3 /2, -H < H < 



(12) 



In Fig. 3 the dependence M(H) is shown for different values of fields H p and H s . It 
is seen that the edge barrier causes peaks to arise on the magnetization curve M(H) at 
H ±H S . With an increase of the ratio H s /H p the amplitude of the peak grows and 
the latter shifts towards higher fields. Such a behavior of magnetization is similar to the 



corresponding dependence reported for wide films in |L5 
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Within the Kim-Anderson model, if H s > H% 2 /H kl , H > H s + H% 2 /H kl , ff > y/2H kl , 
\h[ \ < |if |, the dependence M(ff) is written in the form 



M(H) = 7 < 



H (b 3 - al)/2 + 3H (a - b )/2 + if, 
bl(H - H)/2 - H a 3 /2 + Hb 3 /2 + H + 3H (a - b )/2, 
b 3 (H - H)/2 - H a 3 /2 + 3^ - H )/2 + 3ff a /2, 
Ha 3 + (H~ H )((b*) 3 - a 3 ) + 3H a ((b*) 2 - a 2 )/2- 
3H 2 2 (l-b 2 2a )/2(H kl + \H\)+m 0: 



Hdf < H < ff , 

He x < H < H d f, 

H\ +) <H< H eX) 



0<H <H 



(+) 



Ha 3 + (H~ H )((b*) 3 - a 3 ) + 3H a ((b*) 2 - a 2 )/2- 

3H 2 k2 (l - b 2 2a )/2(H kl + \H\)+m + Hb\j2 + H(l - 3b 2a /2), H^ < H < 0, 

H (bl - al)/2 + 3H (a -h)/2 + H, ff« <H< h[~\ 

b 3 (H -H)/2- H al/2 + H - 3(H + H s )/2 + Hb 3 /2, -H < H < if H, 

(13) 

where m = 6/{ W) j£ 2a yj {y)dy. 

When H < y/2H k i, the dependence M(if) will be fully identical to (11) with the 
appropriate values of parameters a, bi,b 2 ,b 5 , H d f, H ex , iff - * 1 obtained using the KA approach. 

Fig. 4 shows dependence (13) for different values of the parameters related to the Kim- 
Anderson depinning model and edge barrier. Unlike in the Bean model, the magnetization 
curve here has a maximum even without an edge barrier. Given an edge barrier, the M(H) 
peak amplitude increases and the peak itself shifts towards the higher fields. 

B. ac - susceptibility 

If an external magnetic field H changes harmonically in time, H(t) = H cos(uit), mag- 
netization M(t) is also a periodical function with a period T = 2ti/uj, whose Fourier com- 
ponents of the magnetization M(H ) determine the ac-susceptibility harmonics 



Y aC = 
An 



27 



TTffo 



M(Q)e mB dQ = X ' n + i X l 



(14) 
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where 9 = oot. 

It follows from the symmetry of the problem (M(H) = —M(—H)) that all of the even 
harmonics are zero. Fig. 5a-5b shows the imaginary and the real parts of the first (i.e. 
fundamental) and third harmonics of ac susceptibility, calculated by formulae (13) both 
with- and without an edge barrier (curves 1 and 2, respectively). One can see that existence 
of an edge barrier produces a quantitative (first harmonic) and a qualitative (change of sign 
in the real part of the 3rd harmonic) effect on the form of the dependence Xn(Ho)- We believe 
that the sign reversal effect (SRE) which was for the first time predicted by the authors in 



TI| is a vivid manifestation of how an edge barrier affects the Xz{Hq) behavior. It should be 



mentioned that a similar behavior is typical for a specific set of higher harmonics, as well: 
(Ho), x'lk-i(Ho) {k — 2, 3). Note also that the SRE is not observed for the harmonics 
x', k+ i(H ), X 'ik + i(Ho) (k = 1,2). 



IV. GEOMETRY EFFECT ON MAGNETIC CHARACTERISTICS OF TYPE-II 

SUPERCONDUCTORS 

As follows from our analysis, the magnetic characteristics of various-width plates (W < A 
and W ^> A ) in a parallel magnetic field and also of both narrow [W < \ 2 /d) and wide 
(W 3> X 2 /d) films in a transverse magnetic field are quite similar. The analogous similarity 
properties of the magnetization curves were pointed out by [|l8j where curves M(H) were 



compared for macroscopic samples of different geometries posessing no edge/surface 2barrier 



In Fig. 6a the magnetization curves are shown for samples without a bulk pinning. 
Curve 1 was obtained for a narrow film (formula (11) with H p = 0), curve 2 corresponds to 
a wide plate in a parallel field (as calculated within the nonlocal model [0), curve 3 is for 
a wide film in a transverse magnetic field ||,. The magnetic field is measured in units of 
H s , essentially the first- vortex entry field which is different in either case: 
1)H S = %/2^W 0,0, 2)H S = %/2n^W\yd 0, f2TJ, 3)H S = $ q /2^ttA^ @. 
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The magnetization is normalized so as to ensure the condition M(H) = H in the Meissner 
state. As is seen from Fig. 6a, the dependences M(H) in the appropriate units are quali- 
tatively alike, therefore, the magnetization harmonics Xn C {Ho) (Fig- 7a- 10a) also feature a 
similarity. (Note that this property is exhibited by susceptibilities with n > 3, as well). 

A second set of figures, 6b-10b, describes the situations in which account is taken only 
of a bulk pinning (J P (H) = const), and an edge barrier is neglected. Curves 1 correspond to 
narrow films (see dependence (11) at H s = H p ), curves 2 are for wide films in the longitudinal 
geometry [ 16"| , curves 3 were obtained for wide films in the transverse geometry Q . The 



magnetic field is measured in units of H* which for narrow films is equal to j p /Co, for curves 
2 - H* = 2njpW/c, for curves 3 - H* = 4j p d/c . Magnetization is measured in units of 
j p W/8c. Note that in all cases considered the curves M(H) never saturate, which is a direct 
consequence of the nonlocality effect. A similarity between the corresponding dependences 
M(H) obtained in different geometries of the problem is more conspicuous here than in 
the case of a zero-bulk-pinning (see above). Besides, it is obvious that the third harmonic 
(both the imaginary and the real parts) changes sign in transition from the edge-barrier to 
bulk-pinning mode. 

Performing simple calculations (for a solvable case of narrow films/plates) under con- 
dition j p <C j s one obtains the following expressions for ^(Hq) and Xsi^o) in the limit 
H ^>H S 



As follows from above expressions in the case of a superconductor without bulk pinning 
XsiHo H s ) > and x!1{Hq ^> H s ) < 0. However, presence of bulk pinning leads to a 
sign reversal of x^O^o) for sufficiently large amplitudes H > H s (j s /j p ) 3 and of Xz{Hq) for 
amplitudes H > H s (j s /j p )ln(j s /j p ). 
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We wish to point out that for large amplitudes of the magnetic field Hq, the numerical 
results related to the case of a narrow film/plate coincide with the analytically obtained 
asymptotics (15 a,b) for ^(i^o) and ~x^(Hq). For an additional control of the calculation 
accuracy we have also employed a twice smaller numerical step; the results practically have 
not changed (the difference did not exceed 3% ). In our opinion, the above arguments fully 
confirm reliability of our calculations. 

The third situation we have examined corresponds to the j p {H) by the Kim- Anderson 
model (ignoring edge barrier). Magnetization curves and ac-susceptibilities are shown in 
Figs 6c- 10c. Here we consider only one case for wide films, specifically, the longitudinal 
geometry, since there are no analytical expressions to describe the dependence M(H) for 



wide films in the transverse geometry. In [23] the results of numerical calculations of the 



M(H) dependence are reported for a wide film. A similar dependence is obtained (given the 
appropriate choice of parameters) also for a wide plate in a parallel field. Curves 1 in Fig. 
6c-10c correspond to a narrow plate (film), curves 2 are for a wide plate. The magnetic field 
is measured in units of H* which is equal to j p (0)/Co for narrow films and to 2nj p (0)W/c 
for wide films. In both the cases magnetization is measured in units of j p (0)W/8c. The 
dependence M(H) corresponding to curve 1 (Fig.6c) is obtained by substitution of the 
expression H s = H\ 2 j(H^\ + \H\) in (13), while M(H) for curve 2 shown in the same 
figure is obtained using the critical state model for a wide plate, in the same way as it was 
done in p|,|l~3f for an infinite cylinder in a parallel magnetic field. Note a practically absolute 
similarity of the curves at H 3> H*; however, at lower H the differences become noticeable. 
Therefore, at large amplitudes H Q the susceptibility harmonics for these two cases almost 
coincide, whereas at small amplitudes they differ in quantity. 

We should emphasize that the KA model also demonstrates sign reversal in the third 
harmonic of ac susceptibility (only in the real part, though) in transition from the edge- 
barrier to bulk-pinning irreversibility mode. As is shown by numerical analysis of narrow 
samples, in a physically interesting case H > H^i this transition occurs when H% 2 / Hki > 
2H S . The latter corresponds to the condition that the depinning current density in a zero 
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magnetic field j p (0) be twice higher than the density of the edge-barrier suppression current. 
Note that such a situation is practically impossible in superconductors with an ideal barrier 
{js = jcL ~ 10 8 A/cm 2 ). Therefore, the edge barrier effect has to be taken into account in 
study of the magnetic characteristics of type II superconductors and of higher harmonics of 
ac susceptibility. 

V. DISCUSSION 

Let us discuss the applicability conditions for a generalized model of critical state. A 
continuum approximation to describe vortex distribution in narrow samples is valid, if W ^> 
£ (here £ is the coherence length). Introduction of the function n(y) implies essentially the 
averaging of vortex distribution on a scale exceeding the intervortex distance ~ rC 1 ! 2 . It is 
easily shown that this averaging procedure is correct in the field range &o/W 2 < H <C $o/£ 2 - 
For example, the applicability condition of this model (£ <C W < X 2 /d) was met near the 
critical temperature (T — > T c ) for a Sn film in the experiment ||24| . Similar conditions can be 
provided also in HTSC films based, for instance, on YBaCuO with parameters A ~ 0.2/xm 
and d ~ 0.1/im at T « 77K. 

The results for magnetic susceptibilities, obtained in this work are valid on condition 
of quasistatical variation of an external field. The slowness of the magnetic field varia- 
tion implies that there is sufficient time for a quasiequilibrium distribution of current and 
vortex density to set in a sample, and the viscous losses related to vortex motion can be 
neglected. As follows from formula (14), in this case is frequency independent. Absence 
of frequency effects on magnetic susceptibilities was reported in refs PJUfl in the frequency 
domain uj/2n = 30 -j- 1000Hz and at amplitudes Hq = 0.1/iT lmT. 

Consider now the temperature dependence of . Note that the similarity described in 
Section 4 is mostly featured by field dependences Xn C (^o)- The temperature dependences 
of Xn{T) will differ appreciably for samples of different geometries and orientation. This 
is predetermined by a different character of the temperature dependence of fields, H S (T) 
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and H P (T) in wide and narrow samples. Nevertheless, it is quite possible to observe a 
sign reversal on the temperature dependence of X3CO due to edge barrier effect. Indeed, 
let H P (T) ~ (1 - T/T c )°, and H s ~ (1 - T/T c )^ and, besides, a > /? (see e.g. 0,0). 
Then, even if at low temperatures H P (T) 3> H S (T) (bulk pinning dominates), at T ~ T c 
H p (T) <C H S (T), i.e., the edge barrier mechanism of irreversibility will dominate. Such an 
interchange of the irreversibility mechanisms results in a sign reversal of x'31 which takes 
place at a certain crossover temperature T*. Note that a similar crossover of magnetization 
was described in 



It should be noted that a change of irreversibility mechanism may be followed by a sign 
reversal not only in x'siHo), as is demonstrated by the KA model (see Fig. 9a,9c), but also in 
Xs(H ) (within the Bean model; see Fig. 9a, 10a-9b, 10b). This brings up a question about 
how sensitive the ^3(^0) behavior is to a bulk-pinning mechanism. We have considered the 



dependence j p (H) which corresponds to a collective mechanism of vortex pinning [26] in the 
form 



. 7r>0\/ Hk H \H\ > Hh. 

j P o \H\ < H k . 

The behavior of M(H) and Xn C (^o), found using Eq.(16) agrees qualitatively with the 
curves shown in Figs. 6c-10c (the KA model). Apparently, these data should also be valid 
qualitatively for an arbitrary function j p {H) monotonically decreasing with H. Therefore, 
we may conclude that the sign reversal effect in the real part of the third susceptibility 
harmonics, following a change of irreversibility mechanism, is a general property of the 
mixed state, which is only slightly affected by specific features of edge barrier 

Let us discuss now an unexpected similarity of the magnetization curves and ac sus- 
ceptibilities. Indeed, in order to calculate the magnetic characteristics of a specific sample 
one should solve a second-order differential equation (see (A7) in the Appendix) (for a bulk 
superconductor case), the integral Bio-Savart equation (A5) (for a wide film in a perpen- 
dicular geometry) and/or the first-order differential equation (A6) or (A8) (for a narrow 
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film/plate). Surprizingly, the magnetization curves obtained on the basis of the solutions 
of these apparently different equations look quite similar (on condition that the same irre- 
versibility mechanism, i.e. an edge barrier or a specific type of the fluxoid pinning by the 
lattice defects, is employed) |27j . Such geometry-independent similarity obviously follows 
from the fact that current densities and vortex distributions in a narrow superconducting 
film/plate (see Fig. 2) and in samples of other geometries are qualitatively similar in corre- 
sponding ranges of the magnetic field. The difference in samples geometry and size results in 
mere quantitative modifications of the characteristic parameters (e. g. a, b, Hdf, H ex etc.; for 
specific examples see Refs. fpL5| — [T7[| ) and in specific dependences j x (y) n(y) corresponding 
to a sample of selected geometry. Thus, taking into account the above results, a narrow 
film/plate may be suggested as a basic model to study qualitative features of the magnetic 
properties of practically used type-II superconductors. 

The geometry-related similarity revealed here seems to be caused by the one-parametric 
description of each irreversibility mechanism by means of introduction of a phenomenological 
current density for the barrier suppression j s or for the depinning j p {H). The microscopic 
basis of such universality still awaits a more detailed study. 



VI. SUMMARY 

In this work the magnetic characteristics of superconducting plates and films of various 
widths in an external quasistatistical magnetic field H = Hq cos(ut) have been calculated. 
The study was carried out within a generalized model of the critical state, which accounts for 
both edge barrier and bulk pinning. To demonstrate the effect of a surface (edge) barrier, 
particular models were considered, describing individual influence of either irreversibility 
mechanism (edge barrier or bulk pinning) on the magnetization and ac-susceptibility of 
type II superconductors. The obtained results have led us to the following conclusions: 
1. The magnetic characteristics of type II superconductors are determined by the type 
of irreversibility mechanism (bulk pinning or surface/edge barrier) rather than geometrical 
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parameters of samples (plates, films) and their orientation relative to an external magnetic 
field. The latter factors bring about only quantitative changes in the dependences M(H), 
Xlf(Ho), which thus prove to be similar (for a specified mechanism of vortex depinning). 2. 
A sign change effect in the real part of the third (7-th, 11-th) harmonics of the magnetic 
susceptibility is predicted to follow a change of irreversibility mechanism. A generalized 
model of the critical state for narrow films has been used to determine the conditions at 
which a crossover of the xTi^o) dependence takes place. 
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VIII. APPENDIX 

Let us consider thin superconducting film in a perpendicular magnetic field (see fig. lb) 
in a mixed state. The relation between current density and vector potential in the London 
model reads 

Here (p is a phase of the order parameter satisfying the equation 

V x Vy?(r) = 27r5(r-r'), (A2) 

where r' = (x', y') is the vortex coordinates. The Ampere law (in a gauge V • A = 0) 
reads 
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AA = (A3) 

By employing a Green function for the three-dimensional Laplace operator A we rewrite 
(A3) in the integral form 

A(r) = Ao(r) +\ J J j y^ dx ' d y' dz 'i ( M ) 

where Ao(r) is the vector potential of the applied magnetic field H = V x Ao; the inte- 
gration in (A4) is performed over entire sample. With the help of Eqs. (Al)-(A4) we derive 
the Maxwell-London equation for the average current density j(y) = d^ 1 J d ^ 2 j x (y, z)dz in 
a thin-film limit d ^ X 

8nX eff dj(y) 2 r l j(y') , x 

J 7 , + ~ / dy=d [H — n(y)0o . (A5) 

cW dy c 7-i y — y 

Here n(y) is the vortex density and distance is scaled in units of W/2. Equation (A5) 
was derived for the first time in Ref. [p8|1 . For the case of a narrow film X e ff/W ^> 1 when 



the self- field of currents is neglected Eq. (A5) reduces to 

— ^ = H-n(y)<h. (A6) 
c dy 

The numerical solution of the equation (A5) shows that the integral term becomes in- 
significant for sufficiently narrow films W < A e //. In the opposite case of wide films, the 
differential term in (A5) can be neglected (everywhere inside the film except for areas near 
edges). In the latter limit Eq. (A5) reduces to the canonical version of the Bio-Savart 



equation that is conventionally studied in a quasi-two-dimensional situation p9|j30| . 

Consider now the mixed state of a superconducting plate in a parallel magnetic field (see 
fig. 1). By taking curl from both parts of the London relation (Al) with account for the 
relation (A2) and Maxwell equations V x H = 4tt}/c and V • B = we obtain the nonlocal 
equation for the distribution of a magnetic field (see [p~6| , p~T[] ) 
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Hl - A 2 ^ = n(y)0 o , (A7) 
dy l 

where H t is a ^-component of a local magnetic field H; = (0, 0, Hi) (all other components 
vanish due to the symmetry of the problem). In the limiting case A — > (A7) reduces to 
the conventional relation Hi = ufo. 

In case of narrow plates X/W 3> 1 the differential term in (A7) is dominant and, 
moreover, Hi ~ H since the self-field of currents can be neglected. Employing relation 
dHi/dy = An j x /c, valid for longitudinal geometry, we obtain the equation for the current 
density j x 

— ^ = H-n(y)4 . (A8) 
c dy 

Thus, the equations describing the current density distribution for narrow films (see 
(A6)) and for narrow plates (see (A8)) are equivalent. The mathematical equivalence of 
these equations follows from the fact that the current-induced self field may be neglected in 
these two extreme cases. 
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Figure captions 

Fig.l 

Problem geometry: a) plate in a parallel magnetic field; b) film in a perpendicular 
magnetic field. 
Fig.2 

Current densities and vortex distributions in a narrow superconducting film/plate at 
different magnetic fields H. a) < H |< H a , b) H s < H |< H , c) H df < H |< H , d) 
H ex <H[< H df , e)0<Hl< H ex , f) if H < H |< 0, g) -H < H [< ifH 
| indicates increasing field and I indicates decreasing field. 

Fig.3 

Magnetization curves for a narrow thin film with an edge barrier and bulk pinning (Bean 
model): 1) j p = j s , 2) j p = O.Qj s , 3) j p = 0.3 j s , 4) j p = 0. 
Fig.4 

Magnetization curves for a narrow thin film with an edge barrier and bulk pinning (Kim- 
Anderson model): 1) j p (0) = j s , 2) j p (0) = 0.6 3s , 3) j p (0) = 0.3 Js . 
Fig.5 

Dependence of magnetic susceptibilities on the amplitude of a quasistatical applied mag- 
netic field (H = H cos(ut)) of a narrow-film superconductor (Kim- Anderson model of bulk 
pinning: j p {0)/j s = 1/5, H kl = H s )). 
Curve 1 - a sample with only bulk pinning. 
Curve 2 - a sample with edge barrier and bulk pinning, 
a) x'i(Ho), b) x'((Ho), c) ^(#o), d) xS(^o). 

Fig.6 

Magnetization curves for superconductors of various widths and geometries with single 
mechanism of irreversibility. 

a) - superconductor with edge barrier, b) superconductor with bulk pinning (Bean model), 

c) - superconductor with bulk pinning (Kim- Anderson model). 

Curve 1 - narrow film (plate), curve 2 - wide plate, curve 3 - wide film. 
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Fig.7 

Dependences of the real part of the first harmonic of magnetic susceptibility on the 
amplitude of an applied magnetic field H . 

a) - superconductor with edge barrier, b) superconductor with bulk pinning (Bean model), 
c) - superconductor with bulk pinning (Kim- Anderson model). 
Curve 1 - narrow film (plate), curve 2 - wide plate, curve 3 - wide film. 
Fig.8 

Dependences of the imaginary part of the first harmonic of magnetic susceptibility on 
the amplitude of an applied magnetic field H . 

a) - superconductor with edge barrier, b) superconductor with bulk pinning (Bean model), 
c) - superconductor with bulk pinning (Kim- Anderson model). 
Curve 1 - narrow film (plate), curve 2 - wide plate, curve 3 - wide film. 
Fig.9 

Dependences of the real part of the third harmonic of magnetic susceptibility on the 
amplitude of an applied magnetic field H . 

a) - superconductor with edge barrier, b) superconductor with bulk pinning (Bean model), 
c) - superconductor with bulk pinning (Kim- Anderson model). 
Curve 1 - narrow film (plate), curve 2 - wide plate, curve 3 - wide film. 
Fig.10 

Dependences of the imaginary part of the third harmonic of magnetic susceptibility on 
the amplitude of an applied magnetic field Hq. 

a) - superconductor with edge barrier, b) superconductor with bulk pinning (Bean model), 

c) - superconductor with bulk pinning (Kim- Anderson model). 

Curve 1 - narrow film (plate), curve 2 - wide plate, curve 3 - wide film. 
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fig. 2 G.M.Maksimova et. al. 
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fig.5cd G.M.Maksimova et. al. 
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